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Abstract 



We consider a stochastic version of the well-known Blotto game, called the 
gladiator game. In this zero-sum allocation game two teams of gladiators en- 
gage in a sequence of one-to-one fights in which the probability of winning is 
a function of the gladiators' strengths. Each team's strategy consist the allo- 
cation of its total strength among its gladiators. We find the Nash equilibria 
and the value of this class of games and show how they depend on the total 
strength of teams and the number of gladiators in each team. To do this, we 
study interesting majorization-type probability inequalities concerning linear 
combinations of Gamma random variables. Similar inequalities have been used 
in models of telecommunications and research and development. 

Keywords and phrases: Allocation game, gladiator game, sum of exponential 
random variables, Nash equilibrium, probability inequalities, unimodal distri- 
bution. 

MSC 2000 subject classification: Primary 60E15, 91A05; secondary 91A60. 

OR/MS subject classification: Primary: games/group decisions-noncooperative; 
secondary: probability-distribution comparisons. 



1 Introduction 



Borel 



(1921) proposed a game, later dubbed Colonel Blotto game by Gross and 



Wagner ( 1950 ). In this game Colonel Blotto and his enemy each have a given (possibly 
unequal) amount of resources, that have to be allocated to n battlefields. The side 
that allocates more resources to field j is the winner in this field and gains a positive 
amount Qj which the other side loses. The war is won by the army that obtains the 
largest total gain. 

The relevance of Borel precursory insight in the theory of games was discussed in 
an issue of Econometrica that contains three papers by Borel, including the translation 
of the 1921 paper (Borel, 1953), two notes by Frechet ( 1953b|[a ) and one by 
Neumann (1953). 



von 



Borel and Ville (1938) proposed a solution to the game when the two enemies 
have an equal amount of resources and there are n = 3 battlefields. The problem was 
then taken up by several authors, including several other famous mathematicians. 



Gross and Wagner (1950), Gross (1950) provided the solution for a generic n, keeping 
the amount of resources equal and the gain in each battlefield constant (oj = aj). 



Blackett (1954, 1958) considered the case where the payoff to Colonel Blotto in each 



battlefields is an increasing function of his resources and a decreasing function of his 
enemy's resources. Bellman ( 1969 ) showed the use of dynamic programming to solve 



the Blotto game. Shubik and Weber (1981) studied a more complex model where 



there exist complementaries among the fields being defended. In this case the total 
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payoff depends on tlie relative value of capturing various configurations of targets. 



Roberson (2006 ) used n-copulas to determine the mixed equilibrium of the game under 



general conditions on the amount of resources for each player. His analysis is based on 



an interesting analogy with the theory of all-pay auctions (see also Weinstein (2005) 



for the equilibrium of the game and Sahuguet and Persico (2006) for the connection 



between all-pay auctions and allocation of electoral promises) . 



Hart (2008) considered a discrete version of the Blotto game, where player A has 
A alabaster marbles and player B has B black marbles. The players are to distribute 
their marbles into K urns. One urn is chosen at random and the player with the 
largest number of marbles in the chosen urn wins the game. In another version of the 
game, called Colonel Lotto game, each player has K urns where she can distribute 
her marbles. Two urns (one for each player) are chosen at random and the urn with 
the larger number of marbles determines the winner. The discrete Colonel Blotto 
game and the Colonel Lotto game have the same value. In a third version, called 
General Lotto game, given a,b > 0, player A chooses a nonnegative integer- valued 
random variable X with expectation E[X] = a and player B chooses a nonnegative 
integer- valued random variable Y with expectation K[Y] = b. The payoff for A is 
F{X > y) — F{X < Y), where X and Y are assumed independent. The value of the 
game and the optimal strategies are determined. 

Other authors who dealt with the Blotto game and its applications include, for 



instance, Tukey ( 1949 


),Sion and Wolfe (1957 


), 


Friedman (195^ 


I ) , Cooper and Restrepo 


(1967 


), 
), 
), 


Penn (1971), 


Heuer 


(2001), Kvasov ( 


200? 


'), Adamo and Matros (2009), Powell 


(2009 


Golman and Page ( i 


!009 ) and many more. We refer to . 


kovenock and Roberson 


(2010 


Chowdhury, Kovenock, and Sheremeta 




2010 ) for some history of the Colonel 



Blotto game and a good list of references. 

In this paper we deal with a stochastic version of the Colonel Blotto game, called 
gladiator game by Kaminsky, Luks, and Nelson (1984). In their model two teams of 



gladiators engage in a sequence of one-to-one fights. Each gladiator has a strength 
parameter. When two gladiators fight, the ratio of their strengths determines the 
odds of winning. The loser dies and the winner retains his strength and is ready for 
a new duel. The team that is wiped out loses. Each team chooses once and for all at 
the beginning of the game the order in which gladiators go to the arena. 

We construct a zero-sum two-team game where each team also has to allocate 
a fixed total strength among its players. The payoff is linear in the probability of 
winning. We find the Nash equilibria and compute the value of the game. The main 
results are: 

(i) the order according to which gladiators fight has no relevance, moreover knowing 
the order chosen by the opponent team does not provide any advantage; 

(ii) the stronger team always splits its strength uniformly among its gladiators, 
whereas the weaker team splits the strength uniformly among a subset of its 
gladiators; 
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(iii) when the two teams have roughly equal total strengths, the optimal strategy for 
the weaker team is to divide its total strength equally among all its members; 



(iv) when the total strength of one team is much larger than that of the other, the 
weaker team should concentrate all the strength on a single member. 



De Schuymer, De Meyer, and De Baets (2006) consider a dice game that has some 



analogies with ours. Both players can choose one of many dice having n faces and 
such that the total number of pips on the faces of the die is a. The two dice are 
tossed and the player with the highest score wins a dollar. 

The model described below for the probability that gladiator i defeats j, is equiv- 
alent, with different parametrization, to the well-known Rasch model in educational 
statistics, (Rasch, 1980), in which the probability of correct response of subject i to 

for a recent mathematical study of 



see 



Lauritzen 2008 



item j is e"'"^-'' /(I -f- e"'~' 
Rasch models). A similar model has been used also in the theory of contests proposed 
by TuUock (1980), as will be described in Section [7j 

Finding the Nash equilibria of the gladiator game involves an analysis of the 
probability of winning. The key step is a result in Kaminsky et al. (1984) that 
translates the calculation of this probability into an inequality involving the sum of 
independent but not necessarily identically distributed exponential random variables. 

The main theorems are demonstrated through interesting and hard probability 
inequalities, whose proofs are of independent interest and turned out to be more 
complicated than expected. Much of the paper consists of these proofs. We rely 
on Szekely and Bakirov (2003) for some of the technical machinery. The problem is 
cast as a minimization problem involving convolutions of exponential variables and is 
solved by perturbation arguments. A key identity, derived using Laplace transforms, 
directs our perturbation arguments to the analysis of the modal location of Gamma 
convolutions. 

Our inequalities are related to majorization type inequalities for probabilities of 
the form F{Q < t), where Q is a linear combinations of Exponential or Gamma 



variables, that appear in 


Bock, Diaconis, Huffer, and Perlman 


(1987) 


Diaconis and 


Perlmar 


^ ( 


1990 


), 


Szekely and Bakirov 




>003 


) and in 


Telatar 


( 


1999 


), 


Jorswieck and 


Boche 


{' 


2003 


), 


Abbe, Huang, and Telatar 


(2011 


) . The motivation in the last three pa- 



pers, and numerous others, is the performance of some wireless systems that depends 
on the coefficients of the linear combination Q. For stochastic comparisons between 



such linear combinations see Yu (2008, 2011) and references therein. 



Linear combinations of exponential variables appear in various other applications. 



For instance Lippman and McCardle (1987) consider a two- firm model in which learn 



ing is stochastic and the research race is divided into a finite number of stages, each 
having an exponential completion date. The invention is discovered at the completion 
of the A^-th stage. If the exponential times for one firm have parameters that can be 
controlled by the firms subject to constraints, then our results apply to the problem 
of best response and equilibrium allocation strategies for such races. 



4 



Finally, it is well known that the first passage time from to of a birth and death 
process on the positive integers is distributed as a linear combination of exponential 
random variables, with coefficients determined by the eigenvalues of the process' 
generator. For a clear statement, a probabilistic proof, and further references see 



Diaconis and Miclo ( 2009 ) . This allows one to consider R&D type races in which one 



can also move backwards, and applies, for example, to the study of queues, where one 
compares the time until different systems reach a given queue size. 

The paper is organized as follows. In Section[2]we describe the model. In Section|3] 
we determine the Nash equilibria and the value of the game for different values of the 
parameters. Section |4] contains the main probability inequalities used to compute the 
equilibria. Section [5] is devoted to the proofs of the main results. Section [6] deals 
with some monotonicity properties, that follow from our main result and have some 
interest per se. Finally Section [7] contains some extensions and open problems. 



2 The model 

We formalize the model described in the Introduction. Two teams of gladia- 
tors fight each other according to the following rules. Team A is an ordered set 
{Ai, . . . , Am} of m gladiators and team B is an ordered set {Bi, . . . , B^} of n glad- 
iators. The numbers m, n and the orders of the gladiators in the two teams are 
exogenously given. At any given time, only two gladiators fight, one for each team. 
At the end of each fight only one gladiator survives. In each team gladiators go to 
fight according to the exogenously given order. First gladiators Ai and Bi fight. The 
winner remains in the arena and fights the following gladiator of the opposing team. 
Assume that for i < m and j < n ai some point, Ai fights Bj. If Ai wins, the follow- 
ing fight will be between Ai and -Bj+i; if Ai loses, the following fight will be between 
Ai+i and Bj. The process goes on until a team is wiped out. The other team is then 
proclaimed the winner. So if at some point, for some i < m, gladiator Ai fights Bn 
and wins, then team A is the winner. Symmetrically if, for some j < n, Am fights Bj 
and loses, then team B is the winner. 

Team A has total strength ca and team B has total strength cb- The values ca 
and Cb are exogenously given. Before fights start the coach of each team decides 
how to allocate the total strength to the gladiators of the team. These decisions 
are simultaneous and cannot be altered during the play. Let a = (ai, . . . , a^) and 
b = [bi, . . . , bn) be the strength vectors of team A and B, respectively. This means 
that in team A gladiator Ai gets strength and in team B gladiator Bj gets strength 
bj. The vectors a, b are nonnegative and such that 

m n 

^ai = Ca, ^bj = Cb, 

i=i j=i 

namely, each coach distributes all the available strength among the gladiators of his 
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team. 

When a gladiator with strength a fights a gladiator with strength 6, the first 
defeats the second with probability 

(2.1) 



a + h 

all fights being independent. When a gladiator wins a fight his strength remains 
unaltered. The rules of the play and its parameters, i.e., the teams A and B and the 
strengths ca^cb, are common knowledge. Call (^^.^(o., 6) the probability that team 
A with strength vector a wins over team B with strength vector b. 
The above model gives rise to the zero-sum two-person game 

^(m, n, ca, cb) = (m, ca), ^{n, cb), Hm,n) (2.2) 

in which team A chooses a E £^ (m, ca) and B chooses h G e^(n, cb), where 

^(m,CA) = ^{a^,...,am)eWl:Y,ai = c^, (2.3) 
^{n,CB) = l^{h,...,bn) eRl : XJ&i = Cij|, (2.4) 

i^m,n = G^,n ' ^. (2.5) 

The payoff of team A is then its probability of victory Gm „(a, b) minus 1/2. We 
subtracted 1/2 to make the game zero-sum. 



As will be shown in Remark 4.2 below, other models with different rules of en- 



gagement for the gladiators give rise to the same zero-sum game. 



3 Main results 



Consider the game ^ defined in (2.2). The action a* is a best response against b 

if 

a* G argmaxifm„(a, b). 
A pair of actions [a*, b*) is a Nash equilibrium of the game ^ if 
Hm,nia, b*) < if„,„(a*, b*) < if„,„(a*, b), for a\\ a e (m, Ca) and b G ^{n, Cb). 
A pair of actions (a*, b*) is a minmax solution of the game ^ if 

max min Hm,n{o-ib) = min max ifm,n(0', b) = -f^m,n(0'*, b*). 
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Since we are dealing with a finite zero-sum game, Nash equihbria and minmax 
solutions coincide (see, e.g., Osborne and Rubinstein, 1994, Proposition 22.2). The 
quantity Hm,n{cL*, b*) is called the value of the game ^. 

The next theorem characterizes the structure of Nash equilibria of the game 
^(m,n,CA,CB). 



Theorem 3.1. Consider the game '^^{m,n,CA,CB) defined in (2.2). Assume that 
ca < Cb- 

(a) There exists an equilibrium strategy profile {a*,b*) of ^ such that for some J C 
{1, . . . , m} we have 

a*=CA/\J\ for ieJ, a* = for i e J", (3.1) 
b* =CB/n for i G {l,...,n}. (3.2) 

Moreover, all pure equilibria are of this form. 

(b) // 

n 

Cb < -Ca, (3.3) 

n — 1 

then J = {1, . . . , m], so that a\ = ■ ■ ■ = a*^ = CA/m and b\ = • ■ • = b*^ = cb/ti. 

(c) // 

St?/ 

CB > ^^fA, (3.4) 

then J = {i}, that is a* = ca for some z G {1, . . . , m} and o,*j = for all j ^ i, 
and b\ = ■ ■ ■ = b*^ = Cb/u. 

(d) Let to = 1.256431 ■ • • be the root of the equation e* = 1 + 2t. Then for fixed m, 
and Ca and cb such that cb > t^CA, the same conclusion as in ^ holds if n is 
sufficiently large. 



Theorem 3.1 shows that if a vector (a*,b*) is an equilibrium, then so is any 
permutation of a* or b* . Moreover the team with the highest total strength always 
divides it equally among its members, whereas the other team divides its strength 
equally among a subset of its members. This subset coincides with the whole team if 
the total strengths of the two teams are similar, and it reduces to one single gladiator 
if the team has a much lower strength than the other team (see Figures 111 |2l and [3]) . 



FIGURES [TJ § AND |3] ABOUT HERE 

For n = 1, i.e., when team B has a single player, equal strength is always team 
A's best strategy. 
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In order to compute the value of the game ^(m, n, c^, cb), we need the regularized 
incomplete beta function 



J(x, a, (3) 



1 



B(a,/3) 



where 



B(a,/3) = / dt 
When a and /3 are integers, then 

a+/3-l 



r-^{i-tY-^ dt, 
r(a)r(/3) 



(3.5) 



r(a + /3) 



/fx, a. 



j=a V / 



xHl — X 



(3.6) 



For properties of incomplete beta functions see, for instance, Olver, Lozier, Boisverl , 



and Clark (2010). 



Theorem 3.2. Consider the game '^^{m,n,CA,CB)- Assume that ca < cb- 
(a) The value of the game is 



rcB 



rcB + ncA 



(3.7) 



where r is the number of positive a* in the vector a* . In particular 



(b) if (3.3) holds, then the value of the game is 

1 ^ / mcB 



mcB + ncA ' 



(3.8) 



(c) if (3.4) holds, then the value of the game is 

Cb 



1 
2 



Cb + ncA ■ 



1, n 



(3.9) 



In general, to compute the value of the game, one only needs to maximize ( |3.7 ) 
over r = 1, . . . , m; any maximizing r gives an optimal strategy for team A. Figure |4] 
shows the value of the game as Cb varies. Different values of Cb imply different 
numbers of positive a* . 

FIGURE a ABOUT HERE 
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4 Probability inequalities 



x) 



We say that X ~ Exp(l) if X has a standard exponential distribution, i.e., P(X > 
= e~^ for X > 0. 

The main theorems of this paper rely on the following result. 



Proposition 4.1. 'Kaminsky et al. (1984)] The probability Gm,n{0',b) of team A 
defeating B is 



\ i=l i=l 

where Xi, . . . , X^-, ^i, • • • , are i.i.d. random variables, with Xi ~ Exp(l). 



(4.1) 



Remark 4.2. The implication of Proposition 4.1 is that two vectors a, a' of strengths 
that are equal up to a permutation produce the same probability of victory, that is. 



the same payoff function (2.5 ). The same holds for two vectors b, h' . Therefore various 



models, with different rules for the order in which gladiators fight, give rise to the 



same game (2.2). This happens, for instance, in a model where the winning gladiator 
does not stay in the arena to fight the following opponent, but, rather, goes to the 
bench at the end of his team's queue, and comes back to fight when his turn comes. 
This happens also when, at the end of each fight, each coach chooses one of the living 
gladiators in his team at random and sends him to fight. Basically, provided the 
allocations of strength in the two teams is decided simultaneously at the beginning 
and is not modified throughout, any rule governing the order of descent of gladiators 



in the arena leads to the same game (2.2). This is true also for nonanticipative rules 



that depend on the history of the battles so far. The key assumption for this is the 
fact that a winning gladiator does not lose (or gain) any strength after a victorious 
battle. This is parallel to the lack-of- memory property in many reliability models, and 



explains why the probability of winning (4.1) involves sums of exponential random 
variables. 



Note that the main result (Theorem 3.1) does not go through if the allocations 



can also be decided dynamically as battles unfold. In this case the resulting game 
is more complicated and optimal allocations may change according to the observed 
history. For instance consider the case where cb is slightly larger than ca- At the 
beginning, suppose team B spreads the strength uniformly across all its players. If 
team B keeps losing some battles, then it may become optimal to spread the strength 
among only a subset of the surviving players. 



The following theorem is the main tool to prove Theorem |3.1 

Theorem 4.3. Let Xi, . . . , Xm and Yi, . . . , Yn, m,n > 1, be i.i.d. random variables 
with Xi ~ Exp(l). For fixed b> 0, let be as in (2.3) and let 



"15 



G arg min P 

a£j2/ (m,m) 



Y,a,X,<bY,Y, 
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Then 

(a) all nonzero values among a*, . . . , aj^ are equal; 

(b) if m > {n — 1)6, then = ■ ■ ■ = aj^ = 1; 

(c) if m < 2{n — 1)6/3, then a* = m for a single i, 1 < i < m, and a* = 0, forj ^ i. 

5 Proofs of the main results 



The long path to the proof of Theorem 3.1 goes through the following steps: first 



we provide a short proof of Proposition 4.1| for the sake of completeness. Then we 



state and prove three lemmas needed for the proof of Theorem 4.3 Then we prove 



Theorem 4.3, and, resorting to it, we finally prove Theorem 3.1 



Proof of Proposition 4- i First note that if X, Y are i.i.d. random variables with 



X ~ Exp(l), then P(aX > bY) = a/{a + b). Therefore, one can see a duel between 
gladiators i and j as a competition in which the probability of winning is the prob- 
ability of living longer, when their lifetimes are a^X and bjY, respectively. At the 
end of a duel, the winner's remaining lifetime is as good as new by the memoryless 
property of exponential random variables, corresponding to the fact that the strength 
of a winner remains unaltered. The teams' total lives are YliLi ca^i and ^j^j^ 
and the probability that team A wins is that it lives longer, which is Gm,n{o-i b), so 
(O) follows. ' □ 



In order to prove Theorem |4.3| we need several preliminary results. Let Gi,G2, Zi, 
be independent with Gi ~ Gamma(uj, 1), Zi ~ Exp(l), for i = 1,2. For = we 
define = with probability 1. 

Lemma 5.1. Given a\^a\, set a\ = bjux and a2 = 0*2 — blu2. Then 
d 

— P(aiGi + aaGs < x) = {a, - a2)^P(ai(Gi + Z{) + 02(^2 + ^2) < x). (5.1) 
Proof. Let 

F{x) = F{aiGi + 02^2 < x) 

H{x) = P(aiG'i + 02^*2 + aiZi + 02^2 < x) 

and let / and h denote the corresponding densities. Let =Sf denote the Laplace 
transform, that is, 

POO 

^(F) = / e-*^F(x) dx. 
Jo 
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Note that (5.1) is equivalent to 



if 



9_ 
d~5 



Fix] 



[a, - a,)^ ( ^^H{x] 



(5.2) 



Using integration by parts we get 

££ f — = t J e"*^' h{x) dx = t E[exp{-t(aiGi + 02^2 + ai^i + 02^2)}]- 

For the left hand side of (5.2) note that we can interchange differentiation and inte- 
gration, and also that 

= ^(F(x))^log^(F(x)). 
Again by integration by parts we have 

if (F(x)) = ^^if{x)) = ^E[exp{-t(aiGi + 02^2)}]. 



It follows that (5.2) is equivalent to 



1 d 



tdS 

Explicitly this becomes 
1 d 



logi^(/(x)) = (ai - 02) t E[exp{-t(aiZi + 02^2)}]. 



td5 



log[(l + ait)-"^(l + a2t)" 



■"21 



(ai - a2)t(l + ait)-\l + a2t)~\ 



(5.3) 



(5.4) 



Using ai = al + 6/ui, and 02 = — 6/u2, (5.4) is verified by a straightforward 
calculation. □ 



A related result to Lemma 5.1, with a similar type of proof, appears in Szekely 



and Bakirov (2003). 



Lemma 5.2. Given a nonnegative vector (a^, . . . , a"^), let 

ai = a[ + 5 / Ui, a2 = a*2 — 5 /u2^ = a* for 3 < i < m. 

Define 

m n 

Q{a,u) = Y,aiG^-bY,yj, (5-5) 

i=i j=i 

where (a, u) := (ai, . . . , am, Mi, ... , Um), Gi, . . . , Gm, Yi, . . . ,Yn are independent ran- 
dom variables with Gi ~ Gamma(Mj,l), for i = l,...,m and Yj ~ Exp(l), for 
j = l,...,n. Let Zi ~ Exp(l), for i = 1,2 be independent of all other variables. 
Then 

cP' 

— P(g(a, u)<x) = (ai - a2)^^{Q{a, u) + a^Z^ + 02^2 < x). (5.6) 
ab ax^ 
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Proof. Set T = Yl^s ^i^i ~ ^ ^J- Then 



— P(Q(a, u) < x\T) = {a, - as)— P(Q(a, u) + a,Z, + a^Z^ < x|T), 



(5.7) 



which is equivalent to (5.1) with a different x. Talcing the expectation in (5.7) over 
T yields (S. □ 



Lemma 5.3. Let X and Y be independent random variables where Y ~ Exp(l) and 
X has a density f{x) such that 

(i) f{x) is continuously differentiable with a bounded derivative on (— oo,c)o), 

(ii) f{x) > for sufficiently small x G {—oo, oo), 

(iii) f{x) is unimodal, i.e., there exists a G (—00,00) such that f'{x) > if x < a 
and f\x) < if x > a. 

For A > 0, denote the density of X + XY by f\{x). Then f\{x) is unimodal and if 
f'x{xo) = then xq is a mode of fx. Moreover, if \> \q> 0, then any mode of fx{x) 
is strictly larger than any mode of fxo (x) . 



Proof. This result is similar to Szekely and Bakirov (2003, Lemma 1). We provide a 



quick proof using variation diminishing properties of sign regular kernels (see Karlin 



1968). First, since the density of Ay is log-concave (a.k.a. strongly unimodal) its 
convolution with the unimodal f{x) is also unimodal, that is, the pdf of X + XY is 
unimodal (see Ibragimov [l956|[KariSl[l968| . 
Differentiating (justified by (i)) yields 



—00 

-x/X 



fix- z)je-'/^ dz 

/'(z)^e(^-")/^ dz 
A 



l(-oo,.)(^)/'(^)e^/^ dz. 



Suppose fxi^o) = 0. Since f'{z) > for z < a, we know from the representation 
above that fx{x) > if x < a, and hence xq > a. The representation also shows 



that the function e^^"^ fx{x) is nonincreasing in x G 



Therefore f'x{x) > if 



X G (a,xo) and /^(x) < if x > xq. It follows that xq is a mode of fx{x). 
For fixed x, the function l(_oo,x)(-2^)/'( 

function of z does not vanish (by (ii)), 
and has at most one sign cha nge from posi tive to negative (by (iii)), and the kernel 



e^' is strictly reverse rule (see |Karlin||l968|). It follows that / l(-oc,x){^)f'{z)e^/^ dz 



has at most one sign change from negative to positive, as a function of A. Thus, if 
for a given x, /L(x) = and A > Aq, then f'x{x) > 0, and the result follows. □ 
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Proof of Theorem \4.3[ Let Q{a) := Q{a, 1^) as in (5.5). Consider minimizing P((5 (a) < 
0) over 



m 



i=l 



Since fl is compact, and P(Q < 0) is continuous in a, the minimum is attained, say, 
at a* G Q. 



Claim 5.4. In any minimizing point a* of¥{Q < 0) the a* 's take at most two distinct 
nonzero values. Moreover, in the case of two distinct nonzero values, the smaller one 
appears only once. 

Proof. Assume the contrary, say < < a2 < a^. We show below in Case [l] that 
more than two distinct values are impossible by showing that al < 03 leads to a 
contradiction. Similarly Case |2] implies the impossibility of repetitions of the smallest 
of two distinct values. Let ai = al + 6, 02 = 03 — 5, Oj = a*, 3 < z < m. Then by 
(5.6 ) we have 



d_ 
86 



q2 

P(Q(a) <x) = (ai - a2)^P(Q(a) + a^Z^ + a^Z^ < x) 



(5.8) 



where Zi and Z2 are i.i.d. random variables with Zi ~ Exp(l), independent of Q. 
We can focus on x = 0. 



Case 1. al < a^. Since 6 = achieves the minimum, both sides of (5.8) with x = 
vanish at 5 = 0. The density function of Q{cl*) + a^Zi is positive everywhere and 
is log-concave and hence unimodal. By Lemma 5.3, S = Q{cl*) + alZi + a2Z2 has a 
mode at zero. Following Case [2] we show that this leads to a contradiction. 



Case 2. al = a^. Then (5.8) gives 



and 



86^ 



P(g(a) < 0) 



hm^^M^^O 
sio 86 



8^ 



2 lim — P(g(a) + aiZi + 02^2 < x) 
5^0 8x^ 



x=0 



A minimum at 5 = entails 

8^ 



8x'^ 



P(Q(a*) + alZi + alZ2 < x) 



> 0, 



x=0 



showing that S = Q{a*) + alZi + 03^2 has a mode that is nonnegative. 
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Thus S has a nonnegative mode in either case. By Lemma 5.3 and < a^, any 
mode of Q{a*) + a\Zi + a*^Z2 is strictly positive, i.e., 

— P(g(a*) + a[Zi + alZ2 < x) > 0. 



x=0 



The latter expression, multiplied by {a\ — Og) is negative. Using (5.8) with in place 
of Og, however, this implies that ¥{Q{a) < 0) strictly decreases under the perturbation 
{a\, Og) — )■ {a\ + (5, Og — 5) for small 5 > 0, which is a contradiction to the minimality 
at 5 = 0. Note that the crux of the proof is in comparing two perturbations. □ 

Claim 5.5. In any minimizing point a* ofF{Q < 0) the a* 's are either all equal, or 
take only two distinct values, in which case one of them is zero. 

Proof. Assume the contrary, and in view of Claim 5.4, suppose we have 

m 

< al < al = ■ ■ ■ = al_^_i, 1 < k < m, a^_|_2 = • • • = = 0, and i 

at, must be of the form 
k+i ^ (1 + 5)m/(/c + 1), 



m. 



1=1 



Then for some 5 G (0, 1/fc), a* 
a* = (1 - k6)m/{k + 1), al = 
We then have 

k + 1 



0. 



Q{a) = (1 - k6)X + (1 + 6)G - XY, A 



b{k + l) 



m m 
with X ~ Exp(l), G ~ Gamma(/c, 1), Y ~ Gamma(n, 1) independently. We show 
that the minimum oi¥{Q{a) < 0) cannot be achieved in the open interval 5 G (0, 1/A;), 
contradicting the assumption that a* is a minimizer. We have 

P(Q(a) < 0) = P (^1 + 5(1 - (A: + 1)5) < 

where B := X/(X + G). Note that B has a Beta(l, k) distribution, Y/{X + G) has 
a scaled F{2n, 2{k + 1)) distribution, and B and F/ (X + G) are independent. Thus 

,n-l 



P(Q(a) < 0) = CiE 



2/ 



where above and below, Cj > denote constants that do not depend on 5, and 
Di{5) > denote functions of 5 G (0, l/k), and both may depend on other constants 
such as A, k, etc. It follows that 

(1 + 5(1 -(A; + 1)5))"-! 



dF{Q{a) < 0) 
86 



-CiE 



:i-{k+i)B) 



-C2 



+ \A + l + fe)"+fc+i 



(A + 1 + 5(1 - (A: + l)5))'^+'=+i 
(1 + 5x)"-i 



dx 



(5.9) 
(5.10) 
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where 



p = p{6) :-- 



[(A + 1 - 6k){y - 1) - 6kXy] y--\y - 1)'-' dy, 

il + 5)iX + l-5k) 
(X + I + S)il - 5ky 



and (5.10) uses the change of variables 

_ {1 + 6x){X + 1 - 6k) 
^ ~ {X + l + 6x){l-6k)' 

Using the closed form integral 

J' [ky + n{y - 1)] y--\y - if-^ dy = p-{p - 1)^ 



we get 



XSiX + l-Sk) 



X + l + 6 



p--\p - ir~'p'{6) + k{l-{X + my'^-'iy - 1)^-^ dy 



X6{X + l-6k) 



X + l + 6 



{Xn - k)g{6) - A(A + l)p"(p - 1)' 



X + l-5k + Xn5 
D2{6) [k{Xn - k)6^ + {X + l){k- 1)6 + (A + l)(A(n - 1) - k - 2 
(An - k)g{6) 



+ 



Specifically 



A + 1 - 5A; + Xn6 ' 



D2{S) 



(5.11) 



A(5p"(p- 1)^ 



(1 + 5)(A + 1 + 5)(A + 1 - 5A; + Xn6) ' 

It is helpful to determine the sign of g{6) for small 6 > and large 6 < 1/k. Let us 
denote the integral in (5.9) by g{S), which has the same sign as g{6) for 6 G (0, 1/k). 
dor e 

m 



A Taylor expansion yields 



+ 



X {x + k) 



k-l 



dx + o{5) 



{X{n - 1) - k - 2)5 
_k l{X + l)'^+>'+^ ' (A + l)"+fc+2 
= C3{X{n-l) -k-2)6 + o{6), as 5 i 0. 

By direct calculation, 

^(l/A;) = C4(A(n-l)-A;-l). 
We distinguish three cases: 

(i) X{n — 1) > k + 2. Then g{S) > and hence g{S) > for sufficiently small S > 0. 
Moreover, by ( [5ll| ), g'{5) > D3{5)g{5), 6 G (0,1/A;). It follows that g{5) > 
for all 6 G (0, 1/k), i.e., F{Q{a) < 0) decreases in 5 G [0,1/k]. The same holds 
in the boundary case X{n — 1) = k + 2. 



15 



11 



111 



k + 1 < A(n-l) < k + 2. Then g{S) < for sufficiently small S > 0, and g{S) > 
for sufficiently large 6 < 1/k. If the minimum of F{Q{a) < 0) is achieved at 
S* e (0, l/k), then g{S*) = > g'{S*), and g{S) has at least one root in (0,5*), 
say (5**, such that g'{S**) > 0. This contradicts (5.11), however, because the 
term in square brackets strictly increases in 6. 

A(?T,— 1) < k+1. Then g{6) < for both sufficiently large 6 < l/k and sufficiently 
small 6 > 0. Suppose g{6*) > for some S* e (0, 1/A;). If An > A; then a 
contradiction results as in Case (ii). Otherwise the term in square brackets in 
(5.11) is no more than 

(A + l){k - l)k-^ + (A + l)(A(n - 1) - A; - 2) < 0. 

Thus any S G (0, l/k) such that g{6) = entails g'{6) < 0. This is impossible 
as g{6) cannot cross zero from above without first doing so from below. Hence 
g{S) < 0, i.e., F{Q{a) < 0) increases in 5 G [0, 1/A;]. The same holds in the 
boundary case \{n — 1) = k + 1. □ 



We now prove the three statements of Theorem 4.3 



(g) This is an immediate consequence of Claim |5.5l 
(0 Let h{k) = P(Q(a) < 0) with 

1 < k < m, and ak+i 



ai 



ak 



m 



0. 



Comparing ¥{Q{a) < 0) in Case (iii) of the proof of Claim (5.5) at 5 = and 
6 = l/k, we see that if m > b{n — 1), i.e.. 



b{k + l){n- 1) 



m 



<k + l, 



then h[k + 1) < h[k), 1 < k < m. Thus h[k) achieves its minimum at k 



m. 



^ Suppose now m < h{n — 1). According to Case (|i]), if h{k + l){n — l)/m > k + 2, 
i.e.. 



A; + 1 > 



m 



(h{n — 1) — m) 



(5.12) 



then h{k + 1) > h{k). In particular, (5.12) holds for all k if m < 2b{n — l)/3, 
which yields h{m) > ■ ■ ■ > h{2) > h{l), i.e., h{k) is minimized at A; = 1. In 
general h{k) is minimized at some k < \m/{b{n — 1) — m) — 1]. □ 

Proof of Theorem 3.1. (|a]) Using Proposition 4.1 and Theorem 4.3| (a)([b|, once all 
the ai are multiplied by a factor ca/'^i, we can prove that there exists a Nash 
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equilibrium that satisfies (3.1) and (3.2), which we denote as {a*,b*). Assume 
(a, b) is another equihbrium. Because the game is zero-sum, we have 



Hm.,n{0',b) > Hm,nia*,b) > H„ 



and 



Hri 



(a,b*) < H„ 



[a 



(a, b) < H„ 

Thus equahties must all hold. Since b* (equal allocation) is the unique optimal 
response to a*, for the equality to hold in Hm,n{o,*, b) > Hm,n{o,*, b*) we must 
have b = b* . Similarly, for the equality to hold in Hm,nio,, b*) < Hm,n{0'*, b*), 



a must be of the form (3.1). Thus all pure equilibria satisfy (3.1) and (3.2). 



(|b]) Theorem 4.3 |b) guarantees that if = ■ ■ ■ = = ca/tti and 6^ = ■ ■ ■ = 6* = 
cb/u, then a* is the unique best response to b* and vice versa. This proves 
that (a*, b*) is a Nash equilibrium of the game. This equilibrium is unique by 
the argument in part (|a]). 



(|c]) Theorem 4.3 |c]) guarantees that if a* = ca for some i E {1, . . . , m} and a* = 
for all j 7^ i, and = ■ ■ ■ = = Cs/n, then a* is a best response to b* 
and Theorem 4.3 |b) guarantees that b* is the unique best response to a*. This 
proves that {a*,b*) is a Nash equilibrium of the game. Again the argument 
used in part (|a]) shows that all Nash equilibria are of this form. 

([d]) Suppose team A allocates its strength equally among r players, and team B 
adopts the optimal strategy of equal allocation among all n players. Then, as 
n — > oo, the winning probability for team A approaches /(r) := F{cAGr > rcs), 
where Gr is a Gamma(r, 1) random variable. Letting (3 := cb/ca, we get 



/(r) - /(r + 1) 



,p r(r + i) 
1 



dx 



(r+l)/3 



r(r + i: 



r(. + i)^ + ^'^-"^^ 



r(r + 1) 



L^o 



y 



1 + - ) e-y^ f3dy-l 



where we have integrated by parts in the second equality and changed the 
variables y = x/ [5 — r m the third. The integral inside the square brackets 
obviously increases in r. Hence /(r) > /(r + 1) implies /(r + 1) > /(r + 2) > 
■ ■ ■ > f{m). Moreover, if /3 = c^/c^ > to then /(I) > /(2) by direct calculation. 
In this case /(r) is maximized at r = 1 and r = 1 is the optimal strategy for 
team A in the large n limit. □ 
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Proof of Theorem 3.2. (|a]) Using Theorem 3.1 a) we know that for some 1 < r < 
and some permutation vr we have ci*^{i) = ■ ■ ■ = a*(-,^.-) = ca/t, a^(r+i) = ■ ■ ■ 
0'-K{m) = 0, and 6* = ■ ■ ■ = 6* = c^/n. Hence 

m n 

a*Xj ~ Gamma(r, r/c^), ^^^j^' ~ Gamma(n, n/c^). 
i=i i=i 



Therefore, (see, e.g. Cook and Nadarajah, 2006, Cook, 2008) 

rcB 



where / is the regularized incomplete beta function defined in (3.5). 
(|b]) By Theorem 3.1 b) in this case r = m. 
(|c]) By Theorem 3.1 c) in this case r = 1. 



(5.13) 



□ 



6 Monotonicity results 

6.1 Monotonicity of the value 

FIGURE [5] ABOUT HERE 



We mention the following consequence of Theorem 3.1 (see Figure [5]) 



Corollary 6.1. In the game '^^{m,n,CA,CB), if the two teams have equal strength 
(i.e., ca = cb), then the value is positive if m > n, namely, the team with more 
players has an advantage over the other team. Moreover, the value of the game is 
increasing in m and decreasing in n. 

Proof. The team with more players always has the option of not using them all. 
Therefore it cannot be worse off than the team with fewer players. However, since 
equal allocation is the unique best response, using them all is strictly better. The 
same argument proves the monotonicity in m and n. Note that directly verifying this 
from the properties of the incomplete beta function appears nontrivial. □ 

FIGURE [6] ABOUT HERE 



Figure [6] shows an interesting implication of Theorem 3^ team B may be at a 
disadvantage even if ca < cb, and this happens if the number n of its gladiators is 
much smaller than the number m of gladiators in A. As the relative difference in 
strength between the two teams increases, it takes a larger number of gladiators to 
compensate for the lower strength. 
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FIGURES [7] AND E] ABOUT HERE 



As Figures It] and Is] show, if condition (3.4) holds, then team A is at a strong 



disadvantage. The disadvantage increases with the total strength cb and the number 
n of gladiators of team B. The number m of gladiators of team A is totally irrelevant, 
since, in equilibrium, the whole strength ca is assigned to only one gladiator. 



6.2 Related probability inequalities 



If Xi, 



Xfn, and Yi, . . . , 1^ are i.i.d. random variables with Xi ~ Exp(l), and 

mX 



X 



m ^ — ^ n ^ — ^ 

i=i j=i 



mX + nY ' 



then Z has a Beta(m, n) distribution. Hence 

p(x <y) = f(z < 



m 



m + n 



m 



m + n 



m, n 



For m > n, by Corollary 6.1, we have 



F{X <Y)< 



(6.1) 



(6.2) 



Since E[Z] = m/{m + n), is equivalent to P(Z < E[Z]) < 1/2, that is, E[Z] < 
Med[Z]. This is a well known mean-median inequality for beta distributions (see 



Groeneveld and Meeden, 1977) 



The inequality (6.2) has the following interesting statistical implication. If two 



statisticians estimate the mean of exponential variables, and use the sample mean as 
their unbiased estimate, then the statistician with the larger sample tends to have 
a larger (unbiased) estimate. If the two of them bet on who has a larger estimate, 
the one with the larger sample tends to win. For normal variables, or any symmetric 
variables, this clearly cannot happen and P(X < Y) = 1/2. 

Suppose now that the two statisticians share the first n variables, that is, for i = 
1, . . . , n we have Xi = Yi, and the remaining variables Xn+i, ■ ■ ■ , X^ are independent 
of the previous ones. Then 



P(X < Y) 




(6.3) 



By (6.2) the last expression in (6.3) is less than 1/2 if and only if m — n > n, that is, 
m > 2n. It equals 1/2 if m = 2n, and it is larger than 1/2 if m < 2n, in which case 
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(6.2) is reversed. Thus in the bet between the statisticians, if most of the variables 
are in common, the odds are against the one with the larger sample, contrary to the 
previous situation. This was noted by Abram Kagan. 

Our main results can be presented in terms of various other distributional inequal- 
ities or monotonicity. Using (3.6) and Corollary 6.1 we obtain further results that 
cannot easily be proved more directly. We say that X ~ Gamma(a,/3) if X has a 
density 

f{x) = -^e-^^x"-\ x>0. 
T{a) 

Corollary 6.2. For m,n integers the following properties hold: 
(a) The function 

m 

m, n 



m + n 

is decreasing in m for fixed n, and increasing in n for fixed m. 

(b) Let T ~ Binom(m + n — 1, m/(m + n)). Then F{T > m) is decreasing in m and 
increasing in n. 

(c) Let S ~ Poisson(m). Then P(S' > m) is decreasing in m. 

(d) Let R ~ Gamma(m, 1). Then F{R < m) is decreasing in m. 
Proof (@) is a restatement of the last part of Corollary |6.1l 

(|b]) follows from ([a]) and (3.6). 

(|c]) follows from (|b]) by letting — t- oo. 

([d]) follows from ([c]) and the identity 

P(5 > m) = — - / e"* dt. □ 

We say that a random variable Q ~ Geom(p) if P((5i = /c) = (1 — p)^p, k = 
0,1,2,.... 

Proposition 6.3. Let Qi, . . . ,Qm be independent random variables such that Qi ~ 
Geom(l/(l + tti)). Define Q = Y,T=iQi- 

(a) We have 

l-G™,„(a,l„) = P(g<n-l), (6.4) 
where a = (oi, . . . , a^) and 1„ denotes the n- dimensional vector of ones. 



(b) IfYl^i ~ then the probability in (6.4) is minimized when all 's are equal. 
In this case Qi are i.i.d. and Q has a negative binomial distribution. 
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(c) IfE[Q] = m, then E[Q] > Med[Q]. 

Proof. (|a]) The relation ( 6.4[ ) can be explained directly: team A loses if all its 



gladiators together defeat at most n — 1 opponents. Gladiator i from team A 
defeats a geometric random number, Qi, of gladiators of strength 1 from team 
B since he fights until he loses, and he loses a fight with probability 1/(1 + aj). 
Thus if Yl^i Qi — then team A defeats at most n — 1 gladiators altogether, 
and loses. 



(|b]) This follows directly from Theorem 4.3 



Note that E,[Q] = Yl^i'^i- Letting n = m, and using (6.4) and part (|b]), we 
conclude that P(Q < n - 1) > 1 - Gm,n{'^m,W = 1/2. We obtain P(Q < 
E[g]) = P(Q <n)> 1/2, and therefore E[Q] > Med[Q]. □ 



7 Comments and extensions 



The probability in (2.1) is a particular example of contest success functior^ The 
following more general class was considered by TuUock (1980) with the purpose of 
studying efficient rent seeking: 



h^{a, h) 



7 > 0. 



(7.i: 



These functions have been studied, axiomatized, and widely used in different fields 
(see, e.g., Skaperdas 1996 Szymanski 2003, Corchon and Dahm 2010 and many 



others). The reader is referred to Corchon (2007), Garfinkel and Skaperdas (2007), 



Konrad (2009) for surveys on this topic. 



In (7.1), when 7 — >■ oo, then 



h^{a, h) — > hoo{a, b) :- 



ii a > b, 
ii a = b, 
ii a < b. 



This case corresponds to a classical Colonel Blotto situation where the stronger glad- 
iator always wins. If the contest success function hoc is used in our game, then any 
equilibrium strategy for the stronger team assigns the whole strength to one single 
gladiator, and, for < c^, team A loses with probability one and value of the game 
is -1/2. 

In ( 7.1[ ), when 7 — > 0, then 



h^{a, b) —7- h^^a, b) :- 



1 ii a > b = 0, 

ii a,b > 0, 
if = a < 6. 



^ Hirshleifer ( 1989 ) calls it technology of conflict 
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When ho is used as a contest success function in our game, then any equihbrium 
strategy assigns positive strength to every gladiator, therefore in each fight either 
gladiator wins with probability 1/2 and the game reduces to one with two teams of m 
and n gladiators respectively, all having equal power. Then, using (4.1), and (5.13) 
we see that the probability that team A wins is equal to 



Gr 



1 



/ ( ^,m,n 



If ai = ■ ■ ■ = Om = 1, then in (6.4) the random variable Q is negative binomial. Hence 
it is easy to see that 



'"-I /I \ "+J 



j=Q 



n + j — I 
j 



and the value of the game is obtained by subtracting 1/2. 

If the extreme cases 7 = and 7 = cxd are easy to analyze, and the case 7 = 1 
required hard calculations, the remaining cases, i.e., 7 ^ {0, 1, 00} look prohibitive in 
our model. They were considered in easier to deal frameworks by some authors. For 



instance, in a context of rent-seeking, when a contest success function of type (7.1) is 



used. Alcalde and Dahm (2010) show that for 7 > 2 the structure of the equilibrium 



is always the same. 



Friedman (1958) and Robson (2005) consider the case 7 = 1 in a static simul- 



taneous battle context similar to the classical Colonel Blotto model and show that 
the equilibrium strategies for both players involve splitting strength evenly across all 
the battlefields. Roberson (2006) considers the case 7 = 00 and shows that the equi- 



librium mixed strategy of the stronger player stochastically assigns positive strength 
to each battlefield, whereas the one of the weaker player gives zero strength to some 
randomly selected battlefields and randomly distributes the strength among the re- 
maining fields. These results bear some analogy with the structure of the equilibrium 
in our game. 



Tang, Shoham, and Lin (2010) consider contest games where the strengths of 



players are exogenously given and coaches simultaneously choose the order of players 
and then players with the corresponding position fight. This model was used by |Arad 
(|2009|. 
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Figure 2: Number of positive a* as a function of cb for ca = 100, n = 20, and various 
m. 
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Figure 3: Number of positive function of cb for ca = 100, m = 40, and 

various n. 
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Figure 4: Value of ^ as a function of cb G [100,200] for c^i = 100, m = 40, and 
various n. 
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Figure 5: Value of ^ as a function of n for m = 20 and {ca = Cb)- 
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Figure 6: Value of as a function of n for m = 20 and different pairs (c^, c^). 
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Figure 7: Value of ^ as a function of n for = 10 and various cb- 
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Figure 8: Value of ^ as a function of cb > 20 for ca = 10 and various n. 
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